Analytical stripe phase solution for the Hubbard model. 
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■ The self-consistent solution for the spin-charge solitonic superstructure in 
quasi-one-dimensional electron system is obtained in the framework of the 
Hubbard model as a function of a hole doping. Effects of interchain interac- 
tions on the ground state are discussed. Results are used for the interpretation 

1 ' of the observed stripe phases in doped antiferromagnets. 
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It is well known |1J that at low enough temperatures quasi-one-dimensional conductors 
may undergo a Peierls- or spin-Peierls transition and develop accordingly either a long-range 
charge- or spin-order. In the case of a discrete lattice model the commensuration effects 
become important 0. Exact solutions of related Hartree-Fock problems led to the picture 
of a solitonic lattice appearing away from half-filling of the bare electron band @J§]. The 
nature of a soliton was determined then by a corresponding order parameter which was either 
lattice deformation, or density of electronic spin. Nevertheless, recently discovered stripe 
phases in doped antiferromagnets (cuprates and nickelates) f| have attracted attention to 
the problem of coupled spin and charge order pararmeters in the electron systems. Numerical 
mean-field calculations |||7j suggest a universality of the spin-charge multi-mode coupling 
phenomenon in repulsive electronic systems of different dimensionalities. On the other hand, 
those calculations are bound to use small clusters which often makes them inconclusive. 

In this paper an exact analytical solution of the Hartree-Fock problem for a one- 
dimensional electron system at and away from half-filling is presented. This solution pro- 
vides a unique possibility to investigate analytically the structure of the periodic spin-charge 
solitonic superlattice. It also demonstrates fundamental importance of the higher order 
commensurability effects, which result in special stability points along the axis of concentra- 
tions of the doped holes. Though there is no long-range order in the purely one-dimensional 
system due to destructive influence of fluctuations, real cuprates are three-dimensional, 
and therefore, the long-range order survives in the ground state. Hence, we believe, that 
one-dimensional mean-field solutions contain universal features of the stripe phase, which 
are stabilized in higher dimensions. We use derived here single-chain analytical solutions 
as building blocks for the stripe phase in quasi two(three)-dimensional system of parallel 
chains. In this way we have found that short-range (nearest neighbour) repulsion between 
doped holes, in combination with effects of magnetic misfit energy between the chains, nat- 
urally leads to formation of either "half-filled" or "fully filled" domain walls in the low- or 
high doping limits respectively. In both cases these walls separate neighbouring antiphased 
antiferromagnetic domains. The "bi-stripe" solution shows up as well. 

The Hubbard Hamiltonian with the hopping integral t and on-site repulsion U (> 0) 
may be written in the form: 

# = *Ei^+tf£(kM#) a ) (!) 

Here an identity: n^hi^ = ^nf — (S? ) 2 is used. Operators h = n-j- + and S z are fermionic 
density and spin (z-component) operators respectively, and a is spin index. Hamiltonian ([!]) 
has convenient form for Hartree-Fock decoupling in the presence of the two order param- 
eters, i.e. electron spin- and charge-densities , m(x) and p(x) respectively. Single-particle 
eigenstates and eigenvalues of the Hamiltonian Eq. Ql|) in the Hartree-Fock approximation 
can be determined from the Bogoliubov-de Gennes equations derived in : 
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where a = U/4t, the Plank constant is taken as unity, and length is measured in the units 
of the lattice (chain) period a. In these units momentum and wavevector are dimensionless, 
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and velocity and energy posses one and the same dimensionality. The wave function is 
defined in terms of the right- left-movers ^ a ± as 



^<t(x) = ^a+(%) exp(mx/2) + a\I / (J _(x) exp(— iirx/2), 



(3) 



where a =|, j, and the Fermi- momentum is pf = Kp/2. In the case of half-filling the average 
number of electrons per site equals p = 1, i.e. pp = tt/2. The averages m(x) and p(x) are 
defined as: 



(n(x)) = p(x), (S z (x)) = m(x) exp(inx) + m*(x) exp(-iirx). 



(4) 



In the case of half-filling: (S z (i)) = (— l) l m(i), where the order parameter m(x) must be 
real. The total energy is equal to 



W 
2t 



^2 e + dx—(m*(x)m(x) — p 2 (x)), 



(5) 



where e = E/2t, /i is a chemical potential. Next, we introduce p and p as 



p(x) — p + p(x), J p(x)dx = 0, 



and pass to a new basis ^i, ^2, according to ||, which in fact is the Bloch-wave basis in 
the periodic charge- density wave potential Up(x) (in what follows we drop spin index a) : 



^±(x) = exp(=pa / p dx')^ i t2 {x 



(6) 



Using this basis we can rewrite Eq.(|]) in the form 



H 



dx 
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where 



A(x) 



-am(x) exp(2ia / p(x)dx). 



(8) 



It is easy to find from the finite band potential theory || all formal single-soliton solutions 
of the eigenvalue problem (^): 

A(x) = e — iko tanh k^x . (9) 
Here eo is the energy of the localized level, counted from the chemical potential, 

(10) 



kn 



A 2 - f 2 



and 2A is the gap in the energy spectrum. Notice, that we consider A and eo (or equivalent 
pair of variables) as two independent variational parameters. The reason is that, according 
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to Eqs. (0) and (H), there are two independent mean fields "hidden" in A(x). Each of them 
should obey a self-consistency equation. First, we derive equation for m(x), using Eq. (|3|): 

m(x) + m*(x) = (^+(x)^*_(x) + c.c.) (11) 

E<f± 

which can also be rewritten in the ^2(2;) representation Definition given in Eq. (§) leads to 
a self-consistency equation for the variable part of the charge mean-field p(x) = (n(x)) — p : 

p(x) = J] + *+(x)* 2 (x)) - p (12) 

Wave functions of the continuum spectrum are as follows: 

_ ±e =F e + fc + ik tanh fc x ifcE 
2y/lJe(e - e ) - k /L 



the upper (lower) sign on the r.h.s. corresponds to the index "1"("2"); energy e and wave 
vector k are related by equation 



e 2 = Al + k 2 



and L is the length of the system (in units of the lattice spacing). 
Simultaneously, wave functions of the localized state with energy eo become: 



For correct summation over the energy levels in all the equations above, periodic boundary 
conditions are imposed on the wave functions of the continuum spectrum: 

*±(x + L) = V±{x) . 

Then, quantization condition follows: 

k 

kL + arctan — = 2irn, (15) 
k 

where n is an integer. Resulting self-consisting charge and spin components of the single- 
soliton state are obtained: 

p(x) = % \v — 1 H — arcsin — — ], (16) 

2 cosh kox 7r A 

e ± ik tanh k Q x , . 2 . e tanhfc x x 

m[x) = exp[±2za(z/ — 1 H — arcsin -r— J( r)J ex Pl*Xj ? (I'J 

a 7r A 2 L 

where i/ is the filling factor of the localized level [y = 0, 1, 2); and \ is an arbitrary phase. 
As is evident from Eq. (|j), at half filling, i.e. in the commensurate case p — 1, the order 
parameter m{x) must be real. In order to fulfil this condition in Eq. ([IT]) (to the lowest 
order in a for the case v — 0, 2, and precisely, for v = 1) one chooses: e = 0. Eqs. (0), dl7|) 
describe structures of the topological spin-charge solitons (kinks), which are either spinless 
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with charge ±e (single electron charge) at v — 0,2, or chargless with spin 1/2 in the case 
v = 1. In all the cases there are two antiphased antiferro magnetic domains in the system. 
The v = 1 soliton is the stationary excited state of the undoped system. The same is true 
in the cases v = {y — 2), but with one electron removed(added) from(to) the system (i.e. 
still zero doping in the thermodynamic limit). Important is that when doping with holes or 
electrons becomes finite, already the ground state of the system possesses periodic "chain" 
of the alternating spinless solitons and antisolitons of the kind v = or v = 2 (for hole- or 
electron doping respectively), which are described above. This solitonic superlattice, as we 
assume, is a one-dimensional analogue of the stripe phase observed in lightly doped cuprates 
and nickelates || . Existence of the superlattice is related to appearance of the central band 
around e = 0, filled with either spinless "holes" ( v = 0) or spinless "electrons" ( v = 2). 

In order to find a structure of the ground state at finite doping we calculate the total 
energy Eq. (|) of the electron system using solutions ( |TE| ) and (|T7|). We obtain for the 
potential energy: 

^£ = [ dx^(m* (x)m(x) - p 2 (x)) = - — - ^(v - 1 + - arcsin ^-) 2 k - ^Lp 2 . 
It J I lot a o 7r A 2 

(18) 

The other part of the total energy reads: 

-^r=l^ e i = (PF6F + A ln— )+w, (19) 

where e 2 F = p 2 F + Aq, and the term w is of the order L°. 

Minimization of the total energy with respect to A in the order oc L gives the usual 
result § : 

A = 2e F exp(-i), A = — . (20) 

A 7T 

Since our approximation makes sense when A < ep, we conclude that parameter a = 
U/At should be much less than 1 (weak interaction limit). In this limit we see that part in 
the potential energy J ap 2 /2 is of the order a and is much less than other terms, which are 
of the order a -1 or a . 

The total energy of the kink state Eq. can be expressed similar to the work |3| : 

= A [(, - -) costf + I sintf] - ^A (, - ^) 2 sinfl, (21) 

It 71 71 O 7T 

where W is the energy of the state without a soliton, cos# = e /A , sin^ = k /A . As we 
have already seen, eo = for the half-filled case, thus leading to 9 = tt/2. 



Minimization of (|2TD with respect to 6 in the case p ^ 1 (it is assumed that \p — 1| ^> 



Aq/vf) leads to the equation: 

2 7r 2 c\ cv it 

n {9 ~ 2 U)[{1 -3^ )sin9 -^ {e ~2 U) C0S9] = °- (22) 

Solutions are the same as for the Peierls model and are independent of a when o<l. The 
only nontrivial solution is a solitonic excited state with: 
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i/ = l, e ~0, g^O, W - W = 2A /tt; (23) 

where q is the charge of a soliton. This chargeless solitonic excitation corresponds to a 
gradual phase-shift by tt/2 of the argument of the sine function in the periodic ground state 
solution given below in Eq. (^fj). Thus, we see that the structure of the ground state of our 
model should be similar to the one of the Peierls model ||. Hence, we can conclude that for 
the finite hole density = \p — 1| charged (y — 0, 2 ) solitons form a periodic structure 
with the spin period I = 2/\p — 1| equal twice the charge period. 

In the lowest order expansion in a <C 1 this structure is described in terms of elliptic 
functions: 

, s 2K(r')A r r „. A . , 1 , E(r) , . 
A*) = J [sn 2 (A x/r,r) - -(1 - -^)], (24) 

m(x) = AoVksn^Aox/Vkjk), (25) 

where su(Aqx/ \fk, k) is the Jakobi elliptic function with the parameter < k < 1 defined by 
2K{k)\fk/ Aq = l/\p—l\, where K{r) and E{r) are complete elliptic integrals of the first and 
second kind respectively, and r = 2^/k/(k + l), r' = y/l — r 2 . Parameter k varies from k = 1 
at p = 1 where sn(A x) = tanh(A x), to k <^ 1 where sn(A x/^/k, k) ~ sin(7r|p — l|x). 
Simultaneously, E(0)/K(0) = 1, and E{l)/K{\) = 0. 

In the limiting case of "overdoping" : \p— 1| ^> A , in which case k <^ 1 and K(k) m tt/2, 
one has: 

A 2 

m{x) ~ — — - — | sin(7r|p — l\x), (26) 



~ 2| A ° 1|3 cos(27r|p-l|a;), (27) 

in accord with the results of obtained in the limit of main harmonics coupling. 

It is known that in any exactly integrable model there is no commensurability effects at 
the commensurate points JTT| : 

\p -l\=m/n, (28) 

where m, n are relatively prime integers. That is the energy and other system parameters 
continuously depend on the filling factor p. But in our case the term —ap 2 /2 in the potential 
energy @ violates exact integrability and leads to a pinning of the spin density wave m(x). 
As a result at any rational point we obtain a decrease in the total energy of the order |TT| 

Sw oc -ae~ n const . (29) 

To summarise we consider briefly the two-dimensional (2D) case using our ID results 
by adding weak interchain interactions. In the lowest order approximation in the interchain 
hopping integral t± the interaction energy is 

SW = - J ^2 / dx(cOs(fi - (fj) - 1) + Wcoulomb, (30) 
<i,3> ° 
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where J ~ t 2 ± / A , <Pi is the phase of a spin-density m(x) on the 2-th chain, Wcouiomb is the 
Coulomb interaction energy between the charged kinks (solitons). We suppose for simplicity 
that the Coulomb interaction decreases rapidly with the distance and take into account only 
charged kinks on the neighbouring sites of the neighbouring chains, Wcouiomb = N p q 2 /ea', 
where N p is the number of pairs of the charged kinks, e is a dielectric constant, and a' is 
an interchain spacing. In the half-filled system, i.e. in the absence of the charged kinks, 
the minimum of fl30| ) is achieved when (pi = (pj for all i, j, and we have 5W = 0. There 
are two possible ways to create a periodic solitonic structure (solitonic superlattice) in the 
doped 2D system. In the first case the charged kinks would reside on every chain, while in 
the second case only every even (odd) chain would be doped. Compare the energies of these 
two different configurations. 

In the first case we have ipi = <pj and 

SWx = W Coulomb = N h q 2 /2ea, (31) 

where N h is the number of kinks (solitons) and the charge q could be deduced from Eqs. flTop 
and (|24]). This is a 2D stripe pattern with filled "domain walls", i.e. having one spinless 
charged kink per period a' perpendicular to the chains direction. In the second case the 
Coulomb repulsion energy is negligible, but there is an increment in the total energy due 
to the first term in (|30|), i.e. due to a magnetic order misfit: (pi ^ (pj. The minimal energy 
configuration could be achieved now in two ways. One possible pattern corresponds to half- 
filled "bi-stripe" pattern, seen experimentally in some doped manganites ||. Namely, the 
kink - anti-kink pairs of the smallest possible size £ would form on each even(odd) chain in 
order to keep the phase-shift of m(x) with respect to the odd(even) chains equal to 2n. The 
odd(even) chains remain antiferromagnetically ordered. Then the energy is: 

5W 2 ~ N h J£, (32) 

where £ is of the order of the kink width. £ monotonically increases from its value vp/Ao ~ 
e l / x at p = 1, to v F /A ~ e 2/A tan7r|p - l|/2 in the limit \p - 1| > e~ 1/A 0. In the case 
5W\ > 8W2 we have that the alternating kink structure is energetically preferable. Notice, 
that function 5W\ monotonically decreases, while function 8W2 monotonically increases as 
the function of \p— 1|, therefore the sign in the above inequality may change at some sufficient 
doping concentration of the holes. 

An alternative, half-filled single-stripe pattern may arise when J£ 3> 2A/n , i.e. when 
is not too small. In this case the minimum of the energy (|30D is achieved at (pi = (pj. 
For this to be true for any couple of the chains, while Wcouiomb ~ being also fulfilled, the 
doped holes should again reside, say, on every even chain in the form of spinless charged 
solitons v = . But simultaneously, an equal number of the chargeless solitons v = 1 (with 
spin ±1/2) must be formed at all the odd chains in order to maintain the condition ipi = ipj. 
This configuration will be stable if: 

W s N h < £-N h , (33) 
2ea 

where W s = 2 A/ V is creation energy of the chargeless kink. Notice that since the charge q 
monotonically decreases from 1 to as function of doping Nh = L\l — p\ (see Eqs. (|iT)D,(|2"TD), 
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therefore the inequality C|33|) will be not satisfied at high doping densities, and half-filled to 
filled stripe transition would be expected in qualitative accord with experiment 0]. 

Instructive discussions with Jan Zaanen and S.A. Brazovski are highly appreciated by 
the authors. 



8 



REFERENCES 



[1] G. Grtiner, Density Waves in Solids (Addison- Wesley, Reading, MA, 1994). 

[2] I.E. Dzyaloshinskii and A.I. Larkin, Sov. Phys. JETP 34, 422 (1972). 

[3] S.A.Brazovskii, JETP 78, 677 (1980); S. A.Brazovskii, S. A. Gordyunin, N. N. Kirova, 

JETP Lett. 31, 457 (1980). q 
[4] W.P Su, J.R. Schrieffer, and A.J. Heeger, Phys. Rev. Lett. 42, 1698 (1979); K. Machida 

and M. Fujita, Phys. Rev. B 30, 5284 (1984). 
[5] S-W. Cheong et al, Phys. Rev. Lett. 67, 1791 (1991); C.H. Chen, S-W. Cheong, and 

A.S. Cooper, ibid. 71, 2461 (1993); J.M. Tranquada et al, ibid. 78, 338 (1997); J.M. 

Tranquada et al, Nature (London) 375, 561 (1995). 
[6] J. Zaanen and O. Gunnarson, Phys. Rev. B 40, 7391 (1989); V.J. Emery, S.A. Kivelson, 

and H -Q. Lin, Phys. Rev. Lett. 64, 475 (1990); C. Castellani, C. Di Castro and M. 

Grilli, ibid. 75, 4650 (1995). 
[7] H.J. Schulz, J. Phys. (Paris) 50, 2833 (1989); Phys. Rev. Lett. 64, 1445 (1990). 
[8] S.A.Brazovskii, S.I.Matveenko, Sov. Phys. JETP 60, 804 (1984). 
[9] S.I. Mukhin, preprint |cond-mat/9811409[ 
[10] S.I.Matveenko, Sov. Phys. JETP 60, 1026 (1984). 
[11] I.E.Dzyaloshinskii, I.M.Krichever, JETP 56, 908 (1982). 



9 



